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Abstract 



| The spectral flow of the low-lying eigenvalues of the improved and unim- 

proved Wilson-Dirac operator is studied on instanton-like configurations and 
on thermalized quenched configurations at various /3-values and lattice sizes. 
We also investigate the space-time localisation and chirality of the correspond- 
ing eigenvectors. 

<D . 

1 Introduction 

The spectrum of the Dirac operator in QCD around zero is of particular physical 
interest, because in the continuum it is directly related to topological concepts and 
to the phenomenon of spontaneous breaking of chiral symmetry. 

The topological charge of a smooth classical gauge field is related by the Atiyah- 
Singer index theorem to the difference between the number of right- and left-handed 
exact zero modes of the massless Dirac operator in that gauge background. After 
performing the ensemble average of the path integral and in the limit of an infinite 
volume, the Banks-Casher formula relates the chiral condensate (ipip) to the spectral 
density for zero eigenvalues of the Dirac operator fl|]. Moreover, in the framework of 
an effective chiral Lagrangian one finds sum rules for moments of the spectral density 
in individual topological sectors, and the derivative of the spectral density at zero 
can be related to the decay constant of the pion 0. 

The situation is substantially more complicated in numerical simulations of QCD 
on finite discrete lattices. Besides the restriction to a finite volume and the difficulties 
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of denning a topological charge, the use of the Wilson action explicitly breaks chiral 
symmetry. 

In numerical simulations the low-lying eigenvalues of the Wilson-Dirac opera- 
tor are also crucial from a practical point of view. In the quenched approximation 
(almost-)zero modes are believed to be the origin of "exceptional" configurations, on 
which the numerical computation of light quark propagators is difficult and expen- 
sive, and which lead to large statistical errors of hadronic quantities at small quark 
masses ||. In simulations with dynamical fermions the eigenvalues of the Wilson- 
Dirac operator play a central role for the performance of the algorithms and their 
scaling behaviour when approaching the chiral limit ||. 

In this work we investigate the appearance and space-time localisation of zero and 
low-lying modes of the hernitian Wilson-Dirac operator 75M. In particular, we are 
interested in the effects of the Sheikholeslami-Wohlert (SW) improvement term, the 
variation of the zero-mode distribution with the gauge coupling, and the connection 
between the eigenmodes and topological properties of the underlying gauge field. 

This contribution is organised as follows: After the discussion of general proper- 
ties of the spectra of M and 75M in the next section, we investigate in sect. 3 the zero 
modes on simple, instanton-like gauge configurations, before and after heating, and 
with and without improvement. In sect. 4, we present numerical data for thermalized 
quenched configurations on 16 4 lattices, and discuss our results in sect. 5. 

2 Properties of the Eigenmodes of 75M 

The (improved) Wilson-Dirac matrix is given by 

M = 1-k(H+ l -c sw F^a^) , (1) 

where H is the usual hopping term, and k is related to the bare quark mass rriQ by 
k = (8 + 2m )~ 1 . The boundary conditions for the fermion fields are anti-periodic 
in time direction and periodic in the space directions. The SW improvement term 
contains the usual clover combination F^ u of the SU(3) gauge links. 

Due to the 75-hermiticity, M* = 75M75, the eigenvalues of M come in complex 
conjugate pairs \i and /x*. In addition, for csw — the spectrum on even lattices is 
also reflection symmetric along the line Re/i = 1, i.e. 

Mr = (1 - Kp)r Ms = (1 + np)s , 

with the eigenvectors related by the stagger transformation s(x) = (— l)^^ r(x). 

The elementary relation, (n* k — Hi) -r^^Ti = , between eigenvalues /ij and corre- 
sponding right eigenvectors r\ of M, suggests to distinguish two classes of eigenmodes: 

• "bulk" modes with /1 ^ R and vanishing chirality (r, 75?") 

• real modes with /i G R and possibly non-zero chirality 
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All eigenvalues lie within the disk {p = 1 + up : \p\ < 8}. On configurations with 
increasing /3-values, the bulk modes depopulate the vicinity of the real axis except 
for five regions corresponding to the different corners of the Brillouin zone where the 
quark propagator has poles. 

One may expect that some of the real modes are related to topological properties 
of the gauge field by an approximate lattice remnant of the index theorem, but 
of course also topologically trivial gauge fields can have real modes (e.g. U±(x) = 
diag(l, —1, —1) on one time slice and all other U^x) = 1). 

In the following we shall investigate eigenmodes of the hermitian operator 



The normalisation is such that the eigenvalues A(/t) of Q satisfy < A 2 < 1, and we 
refer to the part of the spectrum closest to zero as the "low-lying" eigenvalues. They 
can be computed in an efficient and reliable way by using an accelerated conjugate 
gradient algorithm ||, which also provides approximate eigenvectors and is viable 
for large lattice sizes. 

Recalling the relation M _1 = Q~ lr ys/ (1 + 8k), we note that all hadronic propa- 
gators can be constructed from Q , which has a simple spectral representation in 
terms of the orthonormal eigenvectors of Q. 

The eigenmodes (eigenvalues and corresponding eigenvectors) of Q(k) have a 
non-trivial K-dependence. The real modes, p — 1 — up G R, of M are equivalent to 
zero-modes of Q at K-values k q = 1/p. In general, the spectral flow Q(k) provides 
indirect information about the complex spectrum of M, and in particular the flow of 
the low-lying eigenvalues A(k) probes the spectrum of M close to the real axis. The 
five regions, where bulk modes of M populate the vicinity of the real axis, correspond 
to minima of the low-lying eigenvalues of Q(k) at K-values roughly of the order of 
0.125, 0.25, ±oo, —0.25, and —0.125. In the special case of a flat gauge field each 
eigenvalue p = 1 — up of M corresponds to an extremal value in the spectral flow of 
A = A ■ (1 + 8k) /k at k = 1 /Rep. 

The chiralities of the eigenvectors of Q 



are related to the spectral flow of Aj(k) by the identity 



i.e. the zeroes of the chirality correspond to extrema in the spectral flow of 75M//C. 
Eq. (JD yields the first order k- dependence of X(k), and from the flow of a zero modes 
close to the position Kq, where A vanishes, one can thus estimate k as 



which is very helpful for an efficient iterative search for possible zero modes. 



Q = 75M/ (1 + 8k) . 



(2) 



Xi = (u»,75«i) 



(3) 
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3 Instanton-like Configurations 



We first investigate to what approximation the index theorem holds for Wilson 
fermions on simple instanton-like configurations. These are set up on even lattice 
sizes in the singular gauge as described in |J, followed by a few cooling sweeps 
to reduce boundary mismatches. We consider configurations with up to two well- 
separated instantons or anti-instantons in the same SU(2) subgroup and a net topo- 
logical charge v — 0, ±1, ±2. 

In the physical ^-region we find n + = v zero modes with chiralities of the same 
sign as u, and n_ = with opposite sign chiralities. On multiple-instanton configu- 
rations their positions are split depending on the size and distance of the instantons. 
In the four unphysical regions one has (n + ,n_) = (0,4z/), (6z/, 0), (0,4z/), and (V, 0), 
respectively. In the following we consider only the physical region 0.1 < k < 0.2. 

The spectral flow of the low-lying bulk modes on cold instanton-like configurations 
is characterised by a value of k x which is close to the free-case value of 1/8 and 
essentially independent of the instanton size p. On the other hand, the position 
k,q and the chirality Xo{ K o) °f the zero mode seem to be sensitive to discretisation 
effects and are strongly p-dependent. The following table shows the values for single 
instantions on a 16 4 lattice: 





c sw = 


csw = 1 


p 


«o Xo(^o) 


«o Xo(^o) 


2 


.135 .66 


.1258 .986 


4 


.127 .90 


.1250 .999 


8 


.126 .99 


.1250 .999 



We note the clear effect of the improvement term (using the tree-level value 
csw — 1) : It significantly rises the chirality x( K o) of the zero mode and moves the 
position k towards lower values and closer to k x . 

When the instanton-like configurations are heated by a moderate number of up- 
date sweeps, the zero-mode remains, but the position k changes and fluctuates for 
different heating trajectories. The heating also moves the positions k x , where the 
chiralities of the lowest bulk modes vanish, to values around K ciit (f3), where (3 is the 
value at which the heating was performed. For the improved operator the variation 
of kq is much smaller, and the central value is lower and significantly closer to k x 
than for csw = 0. 

Also on the heated configurations the eigenvectors are centered on the instantons 
(except for k k x ), but the space-time extension of the eigenvector density at k q 
is not as claerly scaling with the size of the instanton as in the cold case. 

4 Quenched Configurations 

To study the spectrum of the Wilson-Dirac operator on quenched configurations, we 
used 10 thermalized 16 4 configurations for each of the gauge couplings, (3 = 5.9, 6.0, 
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Figure 1: Distribution of the zero 
modes on 10 configurations at (3 = 6.0. 
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Figure 2: Typical flow of the chirality 
of the lowest modes. 



and 6.1. The spectral flows were computed for csw — and for the non-perturbative 
value csw {13) of ref. @. 

The distribution of the zero modes for the improved operator at (3 — 6.0 is 
shown in Fig. 1. The peak of the distribution around K cr i t is more pronounced in 
the improved case than for csw = 0. The typical behaviour of the chiralities on one 
of these configurations is illustrated in Fig. 2. The circles denote the lowest (zero) 
modes and the lines connect points from the same bulk modes to guide the eye. One 
notices that the chiralities of the low-lying bulk modes vanish around the same value 

The space-time localisation of the low-lying eigenmodes for the unimproved case, 
is known to change around K cr - lt from a typically exponential decay of the eigenvector 
density to an even stronger, almost point-like, localisation |J. For the improved 
operator we find in the region k <C k x a somewhat stronger localisation of the 
lowest-lying modes than for c$w — 0, but this difference tends to disappear towards 
k x . Above k x a sharp drop of the "participation ratio" indicates again a very strong 
(almost point-like) localisation similar to the unimproved case. The ^-dependence 
of the localisation strength is the same for the lowest-lying bulk modes and for the 
flow of zero modes. Thus, zero modes appearing at very low k < k x are in general 
weaker localised than those at higher k . 

To investigate the localisation of the zero modes in relation to topological prop- 
erties of the underlying gauge field, we approximate the charge density of the cooled 
gauge configurations by a superposition of single-instanton charge distributions 0. 
We find that the peaks closest to the maximum of the eigenvector density are in 
average significantly higher (i.e. correspond to smaller instantons) than the average 
value of the configuration. There also seems to exist a weak correlation between the 
chirality of the eigenmode and the sign of the topological charge of the closest peak. 
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5 Discussion 



On simple instanton-like configurations we confirmed that the improved and unim- 
proved Wilson-Dirac operator approximately satisfy an index theorem similar to the 



case of staggered fermions or QED 2 [10]. The presence of the zero modes is rather 
stable under the effect of roughening the gauge fields by heating. 

Inclusion of the SW improvement term significantly increases the chirality of the 
zero modes and moves k,q closer to « cr it- On thermalized quenched configurations the 
improvement term renders the zero mode distribution stronger peaked around n crit . 
This is in accordance with the picture that in the continuum limit the K-region with 
vanishing mass gap (i.e. with zero modes of Q) shrinks to a point ITT . 



The number of zero modes found on the 16 4 lattices clearly decreases with in- 
creasing (3- values. This, and the much lower abundance of zero modes on 8 4 lattices, 
indicates that the number of zero modes grows with the physical volume. 

The eigenvectors of the improved operator are somewhat stronger localised than 
in the unimproved case, and their positions tend to lie close to peaks of the topological 
charge density of the (cooled) gauge field. The width of these peaks is typically 
smaller than the average for the configuration. 

On the four exceptional configurations encountered by UKQCD at /3 — 6.0 on 
16 3 48 and 32 3 64 lattices we verified the presence of a zero mode at a near K-value 
Ko < K CT it- In all cases the eigenvector is localised close to or on a narrow peak in the 
topological charge density, and the chirality has the same sign as the charge of the 
peak, but is not correlated with the overall topological charge of the configuration. 

The spectral flow of Q(k) can also provide instructive information about the 
complex spectrum of M. For instance, the position k x of the chirality zeroes of 
the lowest bulk modes approximately coincides with K crit (defined as an ensemble 
quantity, e.g. from the pion mass), and may be an interesting quantity to characterise 
the "critical" K-value for individual gauge configurations. 

Further details and numerical data shall be presented elsewhere JT2 . 
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